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Abstract: The use of quantum correlations between photons to separate 
measure even- and odd-order components of polarization mode dispersion 
(PMD) and chromatic dispersion in discrete optical elements is investigated. 
Two types of apparatus are discussed which use coincidence counting of en- 
tangled photon pairs to allow sub-femtosecond resolution for measurement 
of both PMD and chromatic dispersion. Group delays can be measured with 
a resolution of order 0. 1 fs, whereas attosecond resolution can be achieved 
for phase delays. 

© 2013 Optical Society of America 

OCIS codes: (260.2030) Dispersion; (120.3180) Interferometry; (270.0270) Quantum Optics. 



References and links 

[I] H. Kogelnik and R. Jopson, "Polarization Mode Dispersion," in Optical Fiber Telecommunications IVB: System 
and Impairments, I. Kaminow, T. Li, eds. (Academic, Press 2002), pp. 725-861. 

[2] D. Andresciani, E. Curti, E. Matera, B. Daino, "Measurement of the group-delay difference between the principal 
states of polarization on a low-birefringence terrestrial fiber cable," Opt. Lett. 12 844-846 (1987). 

[3] B. Costa, D. Mazzoni, M. Puleo, E. Vezzoni, "Phase shift technique for the measurement of chi'omatic dispersion 
in optical fibers using LED's," IEEE J. Quantum Elect. 18, 1509-1515 (1982). 

[4] CD. Poole, C.R. Giles, "Polarization-dependent pulse compression and broadening due to polarization disper- 
sion in dispersion-shifted fiber," Opt. Lett. 13, 155-157 (1987). 

[5] CD. Poole, "Measurement of polarization-mode dispersion in single-mode fibers with random mode coupling," 
Opt. Lett. 14 523-525 (1989). 

[6] D. Derickson, Fiber Optic Test and Measurement (Prentice Hall, 1998). 

[7] B. Bakhshi, J. Hansryd, PA. Andrekson, J. Brentel, E. Kolltveit, B.K. Olsson, M. Karlsson, "Measurement of 
the differential group delay in installed optical fibers using polarization multiplexed solitons," IEEE Phot. Tech. 
Lett. 11, 593-595 (1999). 

[8] S. Diddams, J. Diels, "Dispersion measurements with white-light interferometry," J. Opt. Soc. Am. B 13, 1120- 
1129(1996). 

[9] P. Williams, "PMD measurement techniques and how to avoid the pitfalls," J. Opt. Fiber. Commun. Rep. 1, 
84-105 (2004). 

[10] D. Branning, A.L. Migdall, A.V. Sergienko, "Simultaneous measurement of group and phase delay between two 
photons," Phys. Rev. A 62, 063808 (2000). 

[II] E. Dauler, G. Jaeger, A. Muller, and A. Migdall, "Tests of a Two-Photon Technique for Measuring Polarization 
Mode Dispersion With Subfemtosecond Precision," J. Res. Natl. Inst. Stand. Technol. 104, 1-10 (1999). 

[12] M.H. Rubin, D.N. Klyshko, Y.H. Shih, A.V. Sergienko, "Theory of two-photon entanglement in type-II optical 
parametiic down-conversion," Phys. Rev. A 50, 5122-5133 (1994). 



[13] D.N. Klyshko, Photons and Nonlinear Optics (Gordon and Breach, 1988). 

[14] O. Minaeva, C. Bonato, B.E.A. Saleh, D.S. Simon, A.V. Sergienko, "Odd- and Even-Order Dispersion Cancella- 
tion in Quantum Interferometiy," Phys. Rev. Lett. 102, 100504 (2009). 



1. Introduction: Dispersion Measurement - Classical versus Quantum 

As optical communication networlcs migrate towards higher 40 Gbps and 100 Gbps data rates, 
system impairments due to dispersion, especially polarization mode dispersion (PMD), become 
a primary issue. This includes not only fiber PMD, but also contributions from switches, ampli- 
fiers, and all other components in the optical path. The fiber PMD and component PMD tend to 
accumulate in different manners as the size of the network grows. In the long length regime, the 
differential group delay (DGD) due to fiber PMD has a known dependance on length, growing 
as VZ HI. In a similar manner, contributions from chromatic dispersion increase linearly in L. 
This known length dependence makes the dispersion of the optical fibers themselves relatively 
straightforward to measure and to take into account. 

In contrast, component PMD was until recently considered to be too small in comparison to 
fiber PMD to affect significant penalties at the system level. Since the introduction of recon- 
figurable add-drop multiplexers (ROADMs), the number of components that could potentially 
contribute to the PMD in a given system has increased significantly. Although the dispersive 
contribution of each separate component is relatively small, together they are capable of accu- 
mulating and of thereby making a significant contribution to the total system impairment. It is 
therefore important to be able to precisely and efficiently measure small values of DGD. How- 
ever, since only fiber PMD was important in the past, no measuring techniques were developed 
for efficient evaluation of small DGD values. With component PMD starting to play a signifi- 
cant role, developing high-resolution evaluation of small PMD values in a single optical switch 
or other small discrete optical component represents a new challenge to optical researchers that 
must be addressed by modem optical metrology. In this paper, we address the measurement of 
dispersive effects in such discrete elements. 

Polarization mode dispersion is the difference between wavenumbers of two orthogonal 
states of light at fixed wavelength, or equivalently, a polarization-dependent variation of a 
material's index of refraction. A number of methods have been developed for measuring it 
EllSllHISllSlTlEim. Many traditional techniques for measuring PMD rely on an interfero- 
metric approach for high-resolution measurements of absolute values of optical delays. This 
approach requires one to use a monochromatic laser source and to keep track of the number 
of interference fringes. Therefore, the accuracy of the approach is limited by the stability of 
the interferometer, by the signal-to-noise level of the detector, and by the wavelength of the 
monochromatic radiation, leading to significant limitations. For example, the use of monochro- 
matic classical polarized light does not allow one to measure the relative delay between two 
orthogonally polarized waves in a single measurement, so several measurements at different 
frequencies must be used to reconstruct the polarization dispersion properties of materials. The 
use of highly monochromatic laser sources creates the additional problem of multiple reflec- 
tions and strong irregular interference that may have detrimental effect on measuring polariza- 
tion dispersion. 

White-light or low-coherence interferometry [|8] is another widely used approach. The ulti- 
mate resolution of such interferometric measurements will depend on the spectral bandwidth 
of the light source. Achieving sub-fs resolution in PMD measurement dictates the use of light 
sources with bandwidth in excess of 200 nm. Generating light of such a bandwidth with a 
smooth spectral profile is not an easy task in itself. Spectral modulations from existing sources 
with bumpy spectra produce 'ghost' features during measurement, leading to complications in 
dispersion evaluation. In addition, the visibility of interference with such super-broadband fight 



is diminished due to dispersion effects. 

Overall, while classical techniques can provide high-resolution measurement of polariza- 
tion mode dispersion they still have limitations in many areas that quantum-based techniques 
can address. For example, entangled photon states intrinsically provide an absolute value for 
polarization optical delay, in contrast to the conventional (classical) case, which is limited to 
determination of delay modulo an integer number of cycles of the light. This is mainly due to 
the fact that quantum interferometry exploits both phase and group velocity effects in the same 
measurement ifTolfTTI . a feat not possible in classical optics. 

The current practical resolution of conventional dispersion evaluation techniques is limited 
to a few femtoseconds (fs). The primary goal here is to use an interferometric setup with an 
entangled photon source to measure the component PMD of a small, discrete optical element 
to sub-femtosecond precision. Ideally, it would be desirable to measure chromatic dispersion 
with the same device, while allowing for the polarization and chromatic effects to be easily 
separable. We will show that this is indeed possible. Due to the frequency-anticorrelation in 
the entangled downconversion source used as illumination, we may independently determine 
the even-order and odd-order parts of the PMD's frequency dependence. Due to the reliance on 
the frequency anticorrelations within pairs of photons, the separation method is intrinsically a 
two-photon quantum effect, and is not present in the classical interferometer 

Classical attempts to simulate this even-odd separation effect by symmetrical chirping and 
anti-chirping of femtosecond laser pulses are constrained to a very narrow size of wavepacket 
thus making it not very practical. The availability of such a separation is useful in a number 
of circumstances. One example is when there is enough pulse broadening (second-order dis- 
persion) to make accurate measurement of group velocity (first order dispersion) difficult. In a 
fiber, group velocity and broadening effects can be separated to some extent by simply taking 
a sufficiently long length of fiber as sample; the longer the fiber, the more accurately each can 
be measured. When dealing with switching elements or other small discrete optical elements, 
this option is not available. Another means must be found to prevent accurate measurement of 
the first-order group delay from being obscured by second-order broadening effects. That is 
what is accomplished here: the location of a dip in the coincidence rate may be used to find 
the group velocity, and this location is unaffected by the second-order broadening as a result 
of the even-order dispersion cancellation. Conversely, although the amount of broadening in 
a single small component may seem negligible, the total broadening from many such compo- 
nents present in a large network may be significant; thus high-accuracy measurements of these 
very small second-order dispersive contributions is important. Separating them off from the 
generally larger first-order contributions makes accurate measurements much easier 

We note that since the component being analyzed is assumed to be relatively small, the prin- 
cipal polarization axes may be assumed to remain constant over the longitudinal length of the 
object and to be independent of frequency, with the dispersive contributions of the two polar- 
ization components remaining independent of each other This greatly simplifies the analysis. 

After a review of background and notation in section |2] three measurement methods will be 
discussed in sections [3]|5] The apparatus of section |3]uses a single detector to make a classical 
measurement; the system is illuminated with a broadband classical light source. In contrast, 
quantum measurements are made using two detectors connected in coincidence with illumina- 
tion provided by a source of entangled photon pairs (spontaneous parametric downconversion, 
(SPDC)). We will examine two quantum measurement setups in sections]?] and |5] In addition, 
in section |5] we give a qualitative analysis that allows the positions of dips (or peaks) indepen- 
dently of the mathematical formalism. 

The two quantum configurations will be distinguished from each other by referring to them 
as type A or type B. They differ only in the presence or absence of a final beam splitter before 



detection, so they may both be implemented in a single apparatus by allowing a beam splitter 
to be switched in or out of the optical path. Similarly, by adding an additional polarizer and 
counting the singles rate at one detector instead of coincidence events, the classical setup may 
also be implemented in the same device. Thus, a single apparatus could be made which is 
capable of performing any of the three types of measurements to be discussed. 

This paper builds on two previous lines of work. The apparatus used for the type A setup was 
introduced previously ifTOl fTTI . where it was shown that quantum interferometry can achieve 
higher resolution than classical methods in measurements of PMD. Separately, the segregation 
of even- and odd-order chromatic dispersion effects was demonstrated in lfT4l . Here, we bring 
the two strands together in a single device (type B), showing that we can separate even- and 
odd-order effects in PMD, as well as in chromatic dispersion, and that we can do so with the 
resolution available to the type A device. 

As a further benefit of the quantum devices over classical methods, note that for the quantum 
cases there is no need to know in advance the principal axis directions of the device or ob- 
ject being measured. Although the incoming photons are aligned along particular axes that are 
linked to a birefringent crystal orientation, their projections onto any rotated pair of orthogonal 
axes (including the principle axes of the sample) will remain equally entangled, allowing the 
method to work without any need to align the axes of the source and the device under test. 

2. Chromatic Dispersion and Polarization Mode Dispersion 



First consider a material for which the index of refraction is independent of polarization. The 
frequency dependence of the wavenumber k = ^^"^^^ is given by a dispersion relation, which 
can be written near some central frequency iio as 

k{Q.o±(o)^ko±a(o + j5a-±Y(o^ + --- (1) 

for |a)| << D.0- The coefficients a, j3, . . . characterize the chromatic dispersion or variation of 
the refractive index with frequency. Explicitly, 
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Rather than looking at the individual terms in the expansion ([T]i, we may also collect together 
all terms containing even powers of (O and all terms containing odd powers to arrive at an 
expansion containing only two terms: 

k{D.o + C0) = keven{,(0)+kacld{(0), (4) 

where 

^,„™(«) = fco + j3«2 + ^(a)4), (5) 

and 

Kdd{(o) = act) + yor' + ff{(o^). (6) 

In the case of nonzero polarization mode dispersion (PMD), the index of refraction varies 
with polarization. We now have two copies of the dispersion relation, one for each independent 
polarization state: 

fe(no±a)) = km±aHO} + ^Hor + ■■■ (7) 



= kH.even{(0)+kH,odd{(0) (8) 

kv{^o±(o) = kv(,±avO} + fivoP' + ■■■ (9) 

kv.even (O)) + kv,c,dd{(0) , (10) 

where H,V denote horizontal and vertical polarization. 

To describe the PMD, we must define quantities that measure the differences between the 
two polarization states; 

Ako^kvo-kno, Aa = av-aH, A/3 = jSy - jS//. (11) 

These parameters are defined per unit length. For the case of primary interest to us, discrete 
fixed-size objects, the formulas should really be written in terms of the relevant lumped quan- 
tities 

A(j) = lAkQ, AA = lAa, AB = /Aj3, (12) 

where / is the axial thickness of the device under study. However, we will continue to use the 
a, j3, and Ako parameters of eq.[TTl both because they are more commonly used, and because 
they allow easy comparison to the formulas used in fiber optics. 

Note that Ako = is a measure of the difference in phase velocity between the two 

polarization modes, while Aa and A)3 are related to the difference in group velocity. Also, it 
should be pointed out that the PMD and the chromatic dispersion are not entirely independent 
effects; in particular, the PMD coefficients themselves {Mq, Aa, A)3) are frequency dependent. 

In the quantum cases, it is convenient to also define T- = DL, where L is the thickness of the 
nonlinear downconversion crystal and D = u^^ — u^^ is the difference of the group velocities 
of the two polarizations inside the crystal. We will restrict ourself to the simplest case of a 
bulk crystal, so the spectral distribution of the downconverted pairs is described by the function 

Glllil 

4>(a)) = sine ^^T-O^ , (13) 

where the sine function is defined by sine {x) = ^^^^ . Photons are emitted from the downconver- 
sion process in frequency- anticorrelated pairs: the frequencies i2o ± ft) in each pair are shifted 
equally, but in opposite directions, from the central frequency Q,o = C0pui„p/2, with the distribu- 
tion of frequency shifts co being given by cS>{(o) of eq. [T3] The downconversion time scale, 
T_, is inversely proportional to the spectral width of the source, and therefore determines the 
precision of the resulting measurements. The spectrum may be made wider by using a thinner 
nonlinear crystal, but this occurs at the expense of reducing the intensity of the downconverted 
light. High intensity and large bandwidth may be obtained simultaneously by use of a chirped 
crystal, although some of the details of the following analysis will then be changed. 

3. Classical PMD Measurement 

An apparatus equivalent to that shown schematically in fig. [T] 18] is commonly used to measure 
polarization mode dispersion. The illumination may be provided by any sufficiently broadband 
light source. For easier comparison with the later sections, we will assume the illumination is 
provided by type II parametric downconversion, but this is not necessary; since we use a single 
detector, the entanglement of the downconverted photons will play no role. 

Assume an arbitrary amount of H and V polarization out of the downconversion crystal, so 
that the incident field in Jones vector notation is proportional to 



(14) 



Half-wave 
plate (45°) 




Non- 



birefringent 
delay 



Polarizer 
at 45° 



Object 
I 



Object n D 
NPBS ^ lJ 



Fig. 1. Classical (single-detector) white-light setup for finding total PMD. 



where A// and Ay are the incoming ampHtudes of the horizontal and vertical components. After 
a horizontal polarizer, we destroy the quantum state and just pick off one component. We can 
think of it as a classical broadband source of horizontally polarized light. 
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For path 1 (lower), the horizontally polarized light accumulates a phase corresponding to the 
path length d\. For path 2 (upper), the horizontally polarized light passes through a ^ wave 
plate with fast axis 45° from the horizontal, converting it into vertically polarized light, 

-' \ f( ^^i<-)\dco=[( Adco. (16) 



1 J J \ J J \ Anico 

In addition, the vertically polarized light in path 2 experiences a phase corresponding to the 
path length d2 and an adjustable delay 8 = cTt- 

At the second beam spUtter, the two components form a superposition of the form 
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Jo = / Ah{co) ( jk{co)id2+S) ) dco, (17) 



with k{(o) = ^ (assuming the paths are in free space). In the absence of any sample after the 
second beam splitter, this superposition will pass through a linear polarizer at 45°, resulting in 

■^0 = 7 -^Jf- y gik{o,)ch j^^ik{co)(d2+S) j d(0, (18) 
The intensity at the detector is then given by 

/:=|y^|2= f \AH{co)\^[l+coii{k{co){M-5))] do). (19) 



Here, hd ~d\~ d2 is the path length difference between the two arms. 

If a birefringent sample of length / is introduced between the last beam splitter and the final 
polarizer, an additional polarization-dependent phase shift is added to the vector in eq. (TT] 

r / gi[k{(0)di+kH(a>)l] \ 

Jo= Ah{(0)\ \ik(o)(d2+S)+kv{m ) do). (20) 




Delay (fs) 

Fig. 2. Interferograms produced by apparatus of fig. [TJ/br samples of different thicknesses. 
For a fixed thickness, the size of the shift may be used as a measure of the difference in 
phase velocities of the two polarizations. 



The resulting intensity at the detector is then: 

/= |7^|2 = y|A//(a))|M+cos(^^(At/-5)-Ayt(to)/) d(0. (21) 

For Type II downconversion, the Ah{(o) and Av(fi)) are both proportional to <t>(a)) = 

sine (jT-to). Plotting eq.|2T|as a function of birefringent delay 5 leads to interferograms such 

as those shown in fig. |2] Each interferogram will be phase shifted (moving the positions of 

the peaks and troughs within the envelope) due to the zeroth order difference in dispersion 

AA:o, while the envelope as a whole will be shifted horizontally due to the first order difference 

in dispersion Aa and broadened due to the second order difference Aj3 . The interferograms 

shown in fig.|2]are shifted by different amounts due to the use of different sample thicknesses. 

In this plot, a 200 nm bandwidth centered at 1550 nm was assumed, with a coherence length of 

„ _ ''•o 
- AA 
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4. Type A Quantum Measurement 

The goal now is to extract the polarization mode dispersion of an object with a higher precision 
than is possible with the classical apparatus of the previous section. In addition, we would like 
to be able to measure the even and odd orders of chromatic dispersion for each polarization. 

The setup ifTOlfTTl is shown in fig. [3] The downconversion is type II so that the two photons 
have opposite polarization {H and V). The photons have frequencies fio ± ft), where 2£Iq is 
the pump frequency. Controllable birefringent time delays Ti and T2 are inserted before and 
after the beam splitter Objects of lengths h and I2 may be placed before and after the beam 
splitter, respectively. Polarizers at angles Q\ and 62 from the horizontal are placed before the 
two detectors. In the following, we will take 0\ = O2 ~ ^ and assume that the beam splitter is 
50/50. Information about which polarization state travels in which branch of the apparatus will 
therefore be erased, allowing interference to occur with maximum visibility. 

Rather than the Jones matrix formalism used in the previous section, it will be more conve- 
nient here to use creation and annihilation operators for horizontally and vertically polarized 
photons. The portion of the output from the downconversion process that is relevant to our 
purposes is the biphoton state 



/fl'ft)<I>(ft))4(£2o + ft))4(no-ft))|0), (22) 



which will serve as the incident state of our setup. The positive-frequency parts of the fields at 
detectors Di and D2, respectively, can be written in the forms 

£{+'(fi) = iy"dto{%(«iy*"('"i)'i+flv(«iy['--^('"i'''+'"i"il}e-'"^ (23) 

f2) - i|^/0){fl,/(t02y'^"('^)(''+'2) (24) 



The coincidence rate is then computed by integrating the correlation function 
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G(2)(fi,f2) 



over the characteristic time scale T of the detectors: 

f-T/l 



(25) 



Rc{Tx,T2)= I dtidt2d^\tut2). (26) 

J-T/2 

Since T is generally much larger than the downconversion time t_, we may safely simplify by 
taking T °o. 

Using eqs. l22l|26l the coincidence rate may be written in the general form ( |[T2l ) 

^.(ti,T2)=/?o{1+CM(ti,T2)}, (27) 

where Rq is a constant (delay-independent) background term and C^^ — J d(0 \<i>{(0)\^ ~ 
The dependence on the time delays is contained in the modulation term 

M(ti,T2) = iy"£/co|<I>(co)|2e-'l^^'('"'-^*^(-'*"l''-2''»^' (28) 

^ |giA<:(-(U)/2+((!^o-fl))T2 _|_g-iAA:((o)/2-((no+fl))T2| 

da)\^{co)\^ (29) 



X 



cos{[Afc(co) -Afc(-a))]/i + Ak{co)l2 + 2coTi + {Q.o + co)t2}, 



where the second form follows by changing the sign of the integration variable in the first term 
of the previous line. It can be seen that even-order PMD terms arising from the pre-beam splitter 
object cancel. Thus, measurements made with the object before the beam splitter will give us 
the odd-order PMD, and measurements made with the object after the beams splitter give the 
total PMD; making both measurements and then taking the difference will provide the even- 
order PMD. We can see the roles of the even and odd parts more clearly by splitting Ak into its 
even and odd parts, then using the identity cos(A +B) = cosAcosB — sinAsinB. The result is: 

A/(ti,T2) = J dco\<$>{co)\^ {cos[Akoddico)i2li + h) + (o{2ti + T2)]cos[Akeyenico)l2 + ^ot:!] 

-sin [AA:„^,/(c(j)(2/i +/2) + fi)(Ti +T2)]sin[Afc„„,(to)/2 + iioT2]}. (30) 

Note that the integrand in the second term is odd in o), so the integral over that term vanishes. 
Therefore, this simplifies to 

M(ti,T2) = yiia)|<I>(a))|^cos[A^„,/,/(a))(2/i+Z2) + t(j(2Ti + T2)] 

X COs[Akeven{CO)l2 + ^0'^2\- (31) 
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Fig. 3. Type A setup for measuring PMD parameters Aa = ay — oih o'jd Aj8 = jSy — j3//. 



We see that the even- and odd-order terms have separated into different cosine terms. 

In the special case that A/3 and all higher order terms vanish, the integral of the previous line 
can be done explicitly: 



2n 

M(ti , T2) = — cos + iioT2] A 



In the last line we have used the result 



Aa(2Zi+/2) + (2Ti + T2) 



where 
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is the unit triangle function. 
The coincidence rate is then 



Rc{t:\ , T2) = /?o i 1 + cos [A^o/2 + i^0T2] A 



Aa(2/i+/2) + (2Ti + T2) 



(35) 



This result is consistent with equation A31 of ilOll . with the caveat that an extra time delay 
Ti has been added here. We now have two possibilities: we can scan over Ti while holding T2 
fixed, or vice-versa. If we scan over X\ with T2 = 0, we find a triangular dip similar to the HOM 
dip, as shown in fig.|4] The first order term in the PMD, Aa shifts the triangular envelope left 
or right, so that the bottom of the dip is at Ti = — ^ (2/i + /2); thus Aa may be determined by 
measuring the location of the minimum. The absolute value of the factor cos (AA:o/2) in front of 
the triangle function gives the visibility of the dip; so measuring the depth of the dip allows Ako 
to be determined. Note that (depending on the sign of cos (AA:o/2)) the "dip" may actually be a 
peak. 

Alternatively, we may scan over T2 while holding X\ — 0. This leads to an oscillating inter- 
ference fringe pattern within the triangular envelope, similar to those of fig. |2] The shift of the 
triangular envelope allows Aa, the first order term in the PMD, to be determined as before. In 
this case, rather than determining visibility, the zeroth order term A^o horizontally shifts the 
fringe pattern by distance T2 ~ within the envelope, allowing determination of Mo from 
the size of this shift. To see clearly the effects of each order of dispersion, fig.|5]shows examples 
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Fig. 4. Scanning over Ti while keeping T2 = 0. The horizontal shift of the minimum away 
from the origin determines AO!, while the depth of the dip determines A^o- The triangle 
function may lead either to a dip (as shown) or to a peak, depending on the sign of the 
cosine. 




Fig. 5. Scanning over Xn while keeping Tj = 0. In (a), a nonzero Aa shifts the envelope 
from its position for Aa = in part (b). The size of the shift can be measured with accuracy 
on the order of 0. 1 fs. In part { c), a nonzero A/tq shifts the locations of the peaks within 
the unshifted envelope. The size of the shift can be measured with accuracy on the order of 
.001 fs = 1 as. 



of such scans in the presence of zeroth-order and first-order dispersion separately. The fringes 
within the envelope as T2 is scanned allow evaluation of the phase delays (the Mq term) to an 
accuracy on the order of attoseconds (10^'^ s) ifTOl . Group delays from the Aa term down to 
the order of 0.1 fs. 

Note that only the differences of the horizontal and vertical polarization parameters (Aa, Aj3, 
etc.) appear in the formulas above. The resulting interferogram is independent of the values of 
the parameters for fixed polarization («//, ay, etc.) and so are insensitive to non-polarization- 
dependent dispersive effects. 

In principle, Fourier transforming experimental data for the coincidence rate and then fit- 
ting parameterized curves to it will allow the determination of higher order PMD parameters. 
However, this requires a large quantity of data to be obtained at high precision. By adding an 
additional beam splitter to the apparatus in the next section, we will arrive at a better method, 
which allows us to extract additional information; namely, it will also give us information about 
the H and V polarizations separately, not just their difference. 

5. Type B Quantum Measurement 

The goal here is to see if additional information may be obtained with a variant of the previous 
apparatus that mixes the final beams via an additional beam splitter This variation is inspired 
by the setup of ref. Ifl4l . in which even and odd portions of the chromatic dispersion were 
separated into different parts of an interferogram, allowing them to be studied independently of 
each other 

Consider the setup in fig.|6] This differs from the setup of the previous section (fig. ID only 
by the addition of an extra beam splitter before the detectors and an additional nonbirefringent 
delay T in one arm, after the first beam splitter Two birefringent samples of lengths l\ and I2 are 
placed before and after the first beam splitter Birefringent delays Ti and T2 are present before 
and after the beam splitter as well, and a nonbirefringent delay T is added to one of the two 
arms after For the sake of definiteness, assume that Ti and T2 delay the vertical (V) polarization 
and leave the horizontal (H) unaffected. The system is illuminated with type II downconversion 
beams. The pump frequency is at 2^0, while the signal and idler frequencies will be written as 
Ho ± CO. We will make use of the fact that the downconversion spectral function is symmetric, 

4>(a)) = 4>(-co). (36) 

We will identify the e and o polarizations with V and H respectively. 

It should be emphasized that in the notation used here, T is an absolute delay, so it must be 
positive. However, Ti and T2 are relative delays of the vertically polarized photon compared to 
the horizontal, and so Ti and T2 may be positive or negative. 

We will find below that the effects of the even and odd orders separate and play different 
roles: the location of each dip in the interferogram (represented mathematically by a triangle 
function in the coincidence rate) is determined by the odd part, while the relative depths of 
the dips are controlled by the even part. We may predict the number and location of each of 
these dips by identifying the ways in which it becomes impossible from the relative timing of 
detection events in the two detectors to identify which photon took which path. To do so, first 
note that the delay between the V and H photons arising before the first beam splitter is 

ATp„ = Ty - T// = Aa /i + Ti . (37) 

(This is the delay due to the object and Ti alone; it is assumed that the intrinsic delay introduced 
by the known birefringence of the crystal itself has been compensated.) There are four possible 
ways in which the delay after the first beam splitter may compensate this pre-beam splitter 




Fig. 6. Type B setup for finding even- and odd-order PMD 



delay, leaving a total delay of zero between the two photons. These are enumerated in the 
table of figure |7] which gives the total post-beam splitter delay AXpost for each case in the final 
column. Setting 

ATpre + ATpo,, = (38) 

for these four possibilities predicts four dips in the coincidence rate at delay values for which 
the difference in the final column vanish; at these values, there is no path information available 
because the two photons arrive at the detector simultaneously, allowing for complete destructive 
interference between paths. 

One additional dip (represented by the second triangle function in eq.|46] below) arises in a 
slightly different fashion. Here the time delay between the two photons is nonzero, but has a 
value that makes identification of the path impossible to identify. When the vertically-polarized 
photon V follows the lower path after the first beam splitter and the horizontally-polarized 
photon H takes the upper, the total relative delay between the two photons is 

Aw,/ = Ty - T// = Aa /i + Ti = Alpre- (39) 

But when the two photons are interchanged (H along the lower path, V along the upper), the 
relative delay is 

Axtotai = Aa [h + h) + (ti + T2). (40) 

If we require these to be negatives of each other (in other words, requiring ls.Zp„st — —2ATpre) as 
in ref. ifTOl . we find the condition Aa {2li 4-/2) + (2Ti + T2) = 0. This leads to AT,o,a/ = —AXpre- 
Because the photons arrive at different times and with different phases we see that in this case 
interference can occur, leading to the possibility that sines or cosines may modulate this term. 
These expectations will be explicitly verified below for the linearized case. 

The coincidence rate may again be written in the general form of eq.|22] The delay-dependent 
modulation term will now be 



M(ti,T2,t) = J £/a)|4>(to)pe-2'''^ie-2iAU,((»);i ^4^^ 

X |l _e-'Pv(^2o+fl))-*:v(no-fl)))/2+2fl)(T+T2)] _^g/p«(no+(»)-<:v(no+(a))/2-(no+(B)T2l 

_^^il(kv(no-m)-kH{na-o>))l2+{^o-o>)z2] _ gi-[fe(no+(B)-<rH(no-<u))/2+2(UT] 
_^i[ikH{Sk,+a>)+kv(£lo-a)))l2+2noT+ino-a>)r2] _ g-il(A://(ao-(B)+/:i/(ao+(H))/2+2noT+(ao+ffl)T2l 
_^^i[{kH{no+a>)-kH[no-a>)-kv[no+a>)+kv{£lo-a>))l2-2a>T2]^^ 
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Fig. 7. TAe post-beam splitter delays corresponding to the four possible outcomes at the first 
beam splitter: one photon can go in each direction, with the vertical following the upper 
(a) or lower (b) path, both photons may follow the lower path (c), or both may follow the 
upper path (d). The second and third columns give the post-BS delays of the vertical and 
horizontal photons, respectively, while the final column gives the difference (vertical delay 
minus horizontal). 



5.1. Linearized dispersion. 

To better understand eq. |4T] let's momentarily ignore the quadratic and higher order terms in 
the dispersion relations. Then we can write: 

ky (fio + «) = kvQ + ay (0 (42) 

kni^Q + co) = kno + aH co (43) 

Ak{Q.o + (o) = Ako+Aa(0 (44) 

Even and odd parts become: 
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^ ay (0 


ky^even 
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This linear approximation allows all of the integrals of eq. |4T]to be done analytically. For 
^(a)) = sine (-^^), repeated use of the integral 

d(a sinc2(fl«)e""(^+^) = -k{ (45) 
a \ 2a J 

leads to the result: 



M(T„T2,T) = -|A(^ j 



(46) 



-4A( ^'^i +Aa(2/i +/2) — ^ ) sin(^ov^2 + iio(T+T2))sin(^o//^2-iioT) 



^ '2 (ti + Aa/i + avh + t+ T2) \ ^ (ti + AaZi - auh - 



2(Ti+Aa(/i+/2) + T2 )^ 

We see that we have five triangle functions centered at the locations predicted in the previous 
subsection; moreover, we see that one of these is modulated by sinesoidal functions as expected. 

A number of specific methods may now be envisioned for extracting the chromatic and po- 
larization mode dispersion parameters from this setup using various combinations of fixed and 
scanned delays. For example, suppose that we scan over Ti, while holding t and T2 fixed. Then 
each of the A factors above gives a triangular spike (of width 2t-) in the coincidence rate cen- 
tered at the value of Ti for which the argument of A vanishes. We can then easily read off the 
locations of these spikes from eq.|46] Explicitly, the various terms of eq.|46]indicate that there 
should be triangular spikes centered at the values 

Ti = -Aa/i, i [tt - Aa(2Zi +/2)] , -{avh + ^ah + ^ + ^2), (47) 
a///2 + T- Aa/i, T2-Aa(Zi + /2)- 

So suppose we have a sample only after the beam splitter {l\ = 0) and then we do three scans 
over Ti, each with different values of t and T2: 

(i) Take | T2 1 large, with t = 0. Then the T2-dependent peaks move far from the origin, off the 
edge of the plot. We will be left with peaks at Ti = and at Ti = OLnh'^ from the locations of 
the latter we can read off Uh- 

(ii) Take |t2| and |t| both large, but with|T2 + t| small. Then we will be left with peaks at 
Ti = and Ti = avh, so we can read off ay. 

(iii) Take |t| large, with T2 = 0. We will be left with peaks at X\ = 0, and Aa/2, so we can 
read off Aa. 

Since it is more difficult to achieve large values for the birefringent delays than for the non- 
birefringent ones, this procedure may not always be the most practical. An alternative version 
will be described below when we examine a special case. 

Finally, notice that some of the triangular spikes will have their heights modulated by cosine 
terms. The arguments of the cosines depend on ^0, so that measuring the heights of these spikes 
relative to the others will allow Ak^ to be determined as well. 

5.2. Adding in quadratic dispersion 

When the quadratic (A)3) term is added back in, analytic expressions can no longer be obtained 
and numerical simulations must be done. An example is shown for one pair of triangular peaks 
in fig. [8] In the figure, unrealistically large values of Aj3 were used to make the effect clearly 
visible. For Aj3 = (red curve), the peaks have the same triangular form predicted earlier As 
Aj3 increased for fixed Aa and NiQ the top of the triangle flattens and gains small oscillatory 
features; the triangle also broadens slightly. 

For realistic values of Ako and Aj3, the alteration of the peak's height by Aj8 is negligible, so 
that the height of the peak can still be used to measure Ak^. The most straightforward method 
to separate the value of Aj3 from A^o is to fit a parameterized curve to the data and look for the 
values of the parameters Ak^ and A)3 that give the best fit. 




Fig. 8. Effect of quadratic dispersion term Aj3 on a pair of triangular peaks. The red curve 
is for Aj8 = 0, the lower ciin'es correspond to increasing values of AP for fixed Aa and 
AkQ. 



Birefringent 
delay 

Non- 
birefringent 
delay 




<S) 



Coincidence 
counter 



Fig. 9. Setup with object and delays only after the first beam splitter. 



5.3. Example: Postponed delay only 

As a special case, we can look at the situation where there is no sample or delay before the 
first beam splitter This is accomplished by setting /j = Ti = 0. (In reality, the downconversion 
crystal itself acts as a birefringent sample before the beam splitter, but a fixed Ti may be inserted 
to cancel it, so that without loss of generality, we can still take the combination Aacnstaih + Ti = 
0.) 

The setup now looks as shown in fig.|9] and the coincidence rate as given by eqs.l27landl46l 

is: 

7? (t, T2) = /?o |2 + 4A J^^^|lj2 ^ ^.j^ j^^^^^ + i2o(T + T2)] sin [koHk - iio^] 
^^ ' 2{avl2 + T + T2) \ f IjaHh + t:)\ , ^ f Aak + ii' 



Holding the remaining birefringent delay T2 fixed and scanning over the nonbirefringent 
delay T, there should be dips at 

T=-avl2 — '^2 and t = -ank, (48) 



as in fig. do] So, running two scans over T using two different values of T2, the location of the 
dip that remains at the same position in both scans gives us the value of an- The other dip 
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Fig. 10. Scan over T for fixed T2> when there is an object only after the first beam splitter. 

moves between the scans; measuring its location during either scan will then give the value of 
ay. Aa is then given by the difference between the two measured values. 

If it is possible to arrange a value of T2 large enough to satisfy T2 = Aa/2, then the term with 
the sines will be large, in which case we may also be able to extract k^n and k^v by scanning 
T over a range of values for which the other triangle functions vanish and fitting the resulting 
data curve to the function sin +iio(T + AaZ2)] sin [^o//^2 — iio^)]. Alternatively, if only 
Nn) is needed (not k^H and kdv separately), it may be simpler to remove the final beam splitter 
(turning the type B apparatus back into type A), then scanning over T2 and find A^o from the 
shift in oscillation fringes via equation eq. [35] 

6. Conclusions 

We have demonstrated that it is possible to separately measure the even and odd-order con- 
tributions to the chromatic dispersion and the polarization mode dispersion, and that such 
measurement have higher accuracy than dispersion measurements accomplished with classi- 
cal white-light interferometry or time-of-flight methods. Summarizing once again the available 
resolutions of the classical and quantum methods, the resolution of group delay measurements 
are limited by (i) the ability to localize the minimum of the envelope in the classical method 
(resolution of a few fs in the plots of fig. 2), and (ii) the downconversion time scale T_ in the 
quantum case (typically tenths of a femtosecond). 

The Type B quantum method provides the high resolution previously available in the type 
A method, but with additional advantage of being able to separate even- and odd-order effects 
for both PMD and chromatic dispersion, without the need for aligning the principal axes of the 
device with the polarization directions of the source photons. In each case, though, separating 
the various orders within the even-order part (separating zeroth order from second order, for 
example) or within the odd-order part is a much more difficult problem, which has yet to be 
solved in a fully satisfying manner As a final point, note that both the type A and type B devices 
are truly quantum, in that a frequency-entangled photon source is required in order for them to 
operate. If a classical source is used, then the frequencies in the two branches will not appear 
in anticorrelated pairs Ho + G) and ii — o (for example in eqs.l28]andl4Tli. so that the required 
cancellations of even or odd orders will not occur in the various terms. 
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